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To study digital Mott insulator LaTiOs and band insulator SrTiOs interfaces, we apply correlated 
band theory (LDA+U) to (n,m) multilayers, 1 < n, m < 9. If the on-site repulsion on Ti is large 
enough to model the magnetic insulating behavior of cubic bulk LaTiOa, the charge imbalance at 
the interface is found in all cases to be accommodated by disproportionation (Ti*+ + Ti^"^), charge 
ordering, and Ti'^^ dxy-orhital ordering, with antiferromagnetic exchange coupling between the spins 
in the interface layer. Lattice relaxation affects the conduction behavior by shifting (slightly but 
importantly) the lower Hubbard band, but the disproportionation and orbital ordering are robust 
against relaxation. 

PACS numbers: 73.20.-r,73.20.Hb,75.70.Cn,71.28.-|-d 



Atomically abrupt ("digital") interfaces (IFs) be- 
tween oxides with strongly differing electronic prop- 
erties (superconducting-ferromagnetic; ferroelectric- 
ferromagnetic) have attracted interest P, due 
to the new behavior that may arise, and for 
Hkely device applications. Hwang and collabora- 
tors 0, 01 have reported coherent superlattices con- 
taining a controllable number of Mott insulator 
[LaTiOa (LTO)] and band insulator [SrTiOa (STO)] 
layers using pulsed laser deposition, with analysis 
suggesting atomically sharp interfaces comparable to 
those produced by molecular beam epitaxy 0] . The 
most provocative result was that the IFs of these in- 
sulators showed metallic conductivity and high mo- 
bility. Electron energy loss spectra (EELS) for Ti 
suggested a superposition of Ti'^"'" and Ti''+ ions 
in the interface region. Incorporating doping with 
magnetic ions, these same materials are being ex- 
plored for spin-dependent transport applications Q . 
Effects of structural imperfections are being stud- 
ied H, B , but the ideal IFs need to be under- 
stood first. 

Single- and three-band Hubbard models with 
screened intersite Coulomb interaction have been ap- 
plied to this IF. Both the Hartree-Fock approxima- 
tion or dynamical mean field with a semiclassical 
treatment of correlation ^3 result in a ferromag- 
netic (FM) metallic IF over a substantial parameter 
range. Ah initio studies reported so far have focused 
on charge profiles [iTl while neglecting correla- 
tion effects beyond the local density approximation 
(LDA) that we address below. Very recent ah initio 
calculations of the effects of lattice relaxation [T^IT^ 



at the IF have provided additional input into the 
Hubbard-modeling of these IFs, allowing the inves- 
tigation of the interplay of correlation effects and 
relaxation in these models. 

The material-specific insight into correlated be- 
havior that can be obtained from first-principles- 
based approaches is still lacking. In LTO/STO su- 
perlattices, the transition metal ions on the per- 
ovskite B-sublattice are identical (Ti) and only 
the charge-controlling A-sublattice cations (Sr, La) 
change across the interface (cf. Fig. QJ. This 
leaves at each IF a Ti02 layer whose local envi- 
ronment is midway between that in LTO and STO. 
In this paper we study mechanisms of charge com- 
pensation at the LTO/STO-IF based on density- 
functional theory calculations [within the general- 
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FIG. 1: Segment of the (1,1) LaTiOa-SrTiOa multi- 
layer, illustrating the cubic perovskite structure (unla- 
beled white spheres denote oxygen). An LaO layer lies 
in the center, bordered by two Ti02 layers, with a SrO 
layer at top and bottom. The lateral size of this figure 
corresponds to the p(2 x 2) cell discussed in the text. 



ized gradient approximation (GGA) "iS]] employ- 
ing the all electron FP-LAPW-method within the 
WIEN2k-implementation including a Hubbard- 
type on-site Coulomb repulsion (LDA-I-U) 0|. We 
focus on (1) the local charge imbalance at the IF and 
its dependence on neighboring layers, (2) the break- 
ing of three-fold degeneracy of the Ti t2g orbitals 
which will be at most singly occupied, (3) magnetic 
ordering and its effect on gap formation, and (4) 
how rapidly the insulators heal (both in charge and 
in magnetic order) to their bulk condition away from 
the IF. 

To explore the formation of possible charge dis- 
proportionated, magnetically ordered, and orbitally 
selective phases at the IF and to probe the relax- 
ation length towards bulk behavior we have inves- 
tigated a variety of (n, m) heterostructures with n 
LTO and m STO layers (1 < n, m < 9), and lateral 
cells of c(2 X 2) or p{2 x 2) The z direction 

is taken perpendicular to the IF. Lattice parameters 
of the systems have been set to the experimental 
lattice constant of STO, 3.92 A, therefore modeling 
coherent IFs on an STO substrate. Bulk STO is a 
semiconductor with a GGA-band gap 15] of 2.0 eV 
(experimental value 3.2 eV), separating filled O 2p 
bands from empty Ti 3d bands. 

Currently LTO (a=3.97 A) and other 3^^ per- 
ovskites are intensively studied because their struc- 
ture is crucial in determining their electronic and 
magnetic behavior 0|. Bulk LTO is an AFM in- 
sulator of G-type (rocksalt spin arrangement) with 
a gadolinium orthoferrite (20 atom) structure; how- 
ever, lattice imaging indicates that only a few layers 
of LTO assume the cubic structure that we use in our 
superlattices. Using the LDA-fU method, an AFM 
insulator is obtained for U > 6 eV, with a mag- 
netic moment Mxi ~ 0.75 /is due to occupation of 
one of the t2g orbitals (orbital ordering arising from 
spontaneous symmetry breaking). FM alignment of 
spins is 50 meV/Ti less favorable. We discuss first 
the (1,1) muhilayer (ILTO/ISTO layer) pictured in 
Fig. n and then consider systems with thicker LTO 
and/or STO slabs to analyze the relaxation towards 
bulk behavior. 

The (1,1) superlattice is modeled in a transverse 
c(2 X 2) cell (not considered in earlier work 13]) with 
two inequivalent Ti ions, which allows disproportion- 
ation within a single Ti layer and is consistent with 
the AFM G-type order in bulk LTO. The on-site re- 
pulsion strength U on Ti was varied from to 8 eV 
to assess both weak and strong interaction limits. 
The Ti moment versus U and the evolution of the 
density of states as a function of U are shown in 
Fig. ^) and b), respectively. Within GGA {U = 0) 
nonmagnetic metallic character is obtained, consis- 
tent with earlier reports. 0, 0] For J7 < 5 eV the 



system is a ferromagnetic metal with equivalent Ti 
ions, i.e. it is qualitatively like earlier results on 
multiband Hubbard models. At J7 w 5.5 eV dis- 
proportionation occurs on the Ti ions, apparently 
weakly first-order as has been found to occur in the 
Naa;Co02 system |l^. Around J7 ~ 6 eV there is 
a half metallic ferromagnet region, but beyond U w 
6.5 eV a gap opens separating the lower Hubbard 
band and resulting in a correlated insulator phase. 
In the following we model the Mott insulating gap 
(0.5 eV) with C/=8 eV. 

The arrangement of disproportionated ions, which 
is charge-ordered (CO) rocksalt, retains inversion 
symmetry and, more importantly, the more highly 
charged cP ions avoid being nearest neighbors. The 
spatial distribution of the occupied d-orbitals in the 
IF Ti02 layer displayed in Fig.Elreveals that besides 
the CO for U > 7 eV this state is orbitally ordered 
(00) with a filled dxy orbital at the Ti'^^ sites, the 
non-degenerate member of the cubic t2g triplet after 
the intrinsic symmetry- lowering effect of the IF. The 
Fermi level lies in a small Mott gap separating the 
occupied narrow d^y band ('lower Hubbard band') 
from the rest of the unoccupied d-orbitals. For ferro- 
magnetic alignment of the spins (Mxi3+ = 0.72^b) 
the gap ensures an integer moment (2.0/iB). 

The system at this level of treatment is a real- 
ization of a quarter-filled extended Hubbard model 







Coulomb Repulsion U (eV) 
2 4 6 8 



_ 0.6 

cn 

3 0.4 

i- 

^ 0.2- 

0.0 
4.0 



S 2.0 
^ 0.0 

LU 



a) 1LT0/1ST0 


HM 


^ ^Ti^: 

FM 

Insulator 



^--'^^FMMetal 



-2.0 





I 








■ 1 1^ 




s 


u=o 


eV 


U=5 


eV 


U=6 


1 

eV 


35 


eV 



DOS (arb. units/eV,spin) 



FIG. 2; a) Phase diagram of the Ti moments for the 
(1,1) superlattice in a transverse c(2 x 2) cell, versus the 
on-site Coulomb repulsion strength U on the Ti 3d or- 
bitals. b) Density of states (spin direction indicated by 
arrows) of the (1,1) superlattice for different U values. 
Disproportionation occurs in a weak first-order manner 
around U «5.5 eV. HM FM indicates a region of half 
metallic ferromagnetism before the Mott gap appears 
around U « 6.5 eV. 
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FIG. 3: 45° checkerboard charge density distribution 
of the occupied 3d states in the charge-ordered Ti02 
layer in the FM (1,1) multilayer. Orbital-ordering due 
to dxy orbital occupation is apparent. The positions of 
0-, Ti^+ and Ti*"'"-ions are marked by white, dark blue 
(black) and light blue (grey) circles, respectively. 



(EHM) system. The Hubbard model itself is metal- 
lic at quarter-filling; when intersite repulsion is in- 
cluded [23, 1^ it becomes CO and insulating. The 
intersite repulsion is included correctly in first prin- 
ciples methods and that combined with the on-site 
repulsion {U) gives charge ordering. 

The calculation was extended to a larger p(2 x 2)- 
cell to allow antiferromagnetic alignment of the 
Ti'^+ spins. We obtain the same CO/00 state with 
an occupied d^jy-orbital on every second IF Ti ion, 
giving a checkerboard ordering of Ti'^^ and Ti''^+, 
regardless of whether the spins are aligned or an- 
tialigned. AFM coupling is preferred by 80 meV per 
p{2 X 2)-cell for the (1,1) superlattice (a spin-spin 
exchange coupling of | J| = 10 meV). For heterostruc- 
tures containing a thicker LTO slab, however, AFM 
coupled spins on the 50% diluted p(2x2) mesh in 
the IF layer will not match the AFM G-type order 
on the LTO side of the slab, where spins in the IF-1 
layer couple antiparallel with a c(2x2)-periodicity. 
Due to this frustration, AFM alignment within the 
IF layer may become less favorable. 

The (n,m) superlattice. To examine charge- and 
spin-order relaxation towards bulk behavior, and ob- 
serve charge accommodation at more isolated IFs, 
we have studied several thicker slabs containing 
{n,m) layers of LTO and STO, respectively, with 
1 < n,TO < 9. Following the experiment of Ohtomo 
et al. we present results specifically for the (1,5) 
and (5,1) as well as the (5,5) superlattices, all of 
which we find to be disproportionated, CO and 00, 
and insulating in the strong interaction regime. As 
is clear both from the layer-resolved magnetic mo- 
ments presented in Table |l] and the layer resolved 
projected DOS for the (1,5) superlattice in Fig.^), 



TABLE I: Layer resolved magnetic moments (in /is) of 
the Ti ions in (n, m) superlattices. Due to the c(2 x 2)- 
lateral unit cell there are two inequivalent Ti-ions in each 
layer, (n, m) denotes a multilayer containing n LTO and 
m STO layers. The IF moments are nearly bulk-like 
and become so at the layer next to the IF layer. (1,5)* 
denotes a configuration where the interlayer distances 
were relaxed according to Ref . IT^ . 



System 
(n, m) 



LTO 



IF-2 



IF-1 



IF 

IF 



STO 

IF-fl IF+2 



(1,1) - - 0.72/0.05 

(1,5) - - 0.71/0.05 0.0/0.0 0.0/0.0 

(1,5)' - - 0.50/0.08 0.0/0.0 0.01/0.01 

(5,1) 0.73/-0.73 -0.73/0.73 0.70/0.05 

(5,5) 0.73/-0.73 -0.73/0.73 0.70/0.06 0.0/0.0 0.0/0.0 



the IF Ti02 layer, and only this layer, is CO/00 
with Ti3+ and Ti^+ distributed in a checkerboard 
manner. At every second Ti-site the t2g states 
split according to the IF-imposed symmetry low- 
ering, and the d^y orbital becomes occupied. The 
t2g states on the Ti'^^ ions remain essentially de- 
generate, and there is only a tiny induced moment 
M'pi4+ = 0.06/xb. Ti ions in neighboring or deeper 
layers in the STO part of the slab have the con- 
figuration 3dP and are nonmagnetic, while those on 
the LTO side of the slab have the configuration 
and are AFM G-type ordered. Thus the charge 
mismatch is localized at the interface layer, with 
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FIG. 4: Layer-resolved density of states of a a) struc- 
turally ideal and b) relaxed (1,5) multilayer. The two 
topmost panels show Ti'^"'" and Ti^"*" at the IF, the suc- 
ceeding panels show the behavior of the Ti-ion in deeper 
layers of the STO part of the slab. While for the ideal 
geometry rapid relaxation of the electronic structure to 
bulk form versus distance from the interface (IF) takes 
place, in the relaxed structure the electronic relaxations 
involve deeper lying layers. 
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bulk LTO and STO character quickly re-emerging 
on neighboring layers. Consequently, these results 
indicate a relaxation length much less that the 1- 
2 nm value estimated from the EELS data The 
same CO/00 results have been obtained on a vari- 
ety of (n, m) LTO-STO superlattices; these repeat- 
edly emerging insulating ordered IF phases are very 
robust. 

However, the systems discussed so far are struc- 
turally perfect with ideal positions of the atoms in 
the perovskite lattice. In the following we discuss 
the influence of lattice relaxations on the electronic 
properties of the system. Recently, two DFT studies 
using GGA and the LDA+U approach 13] inves- 
tigated structural relaxations in LaTiOa/SrTiOs su- 
perlattices, finding that Ti-ions at the IF are dis- 
placed by O.I5A with respect to the oxygen ions lead- 
ing to a longer Ti-Ti distance through the LaO layer 
than through the SrO-layer. This "ferroelectric" -like 
distortion decays quickly in deeper Wing layers. Us- 
ing the relaxations reported in Ref . [13 , we repeated 
the calculations for the (l,5)-heterostructure. The 
resulting layer-resolved projected DOS at the Ti- 
ions is displayed in Fig. ^). The most prominent 
feature is that for the relaxed structure the dxy-ha,nd 
(the lower Hubbard band) has been shifted up by 0.4 
eV, leaving it incompletely (~70%) occupied. The 
charge is distributed in the minority spin channel at 
the Ti'^^-sites (hybridization with 02p bands) re- 
ducing the magnetic moment from 0.71/^b in the 
ideal structure to 0.50//B- Additionally there is a 
small contribution to conductivity of Ti*+ in deeper 
lying layers in the SrTiOa-host whose d-bands now 
slightly overlap the Fermi level. Hence it is the lat- 
tice relaxations that result in a metallic heterostruc- 
ture and a longer healing length towards bulk be- 
havior, in agreement with the experimental observa- 
tions 3, 4] in spite of a majority of the charge being 
tied up at the IF. Still, the CO/00 arrangement 



remains; it is robust with respect to relaxation and 
tetragonal distortion. 

Now we summarize. While the behavior of the 
many superlattices that we have studied produce ro- 
bust results for the IFs which is easily understood, 
they have a charge- and orbital-ordered character 
that was unanticipated from the original reports on 
these heterostructures. If the interaction strength 
U within the Ti 3d states is large enough to repro- 
duce the AFM insulating state in cubic LTO, then it 
is more advantageous for the local charge imbalance 
to be accommodated within the IF layer itself, which 
can be accomplished by disproportionation, followed 
by charge order with Ti'^+ and Ti*+ distributed 
in a checkerboard manner. The interface layer is 
orbitally-ordered, with an occupied d^jiy-orbital; this 
symmetry breaking is due to the intrinsic IF symme- 
try and is unaffected by atomic relaxation. Indeed 
both disproportionation and orbital ordering are in- 
sensitive to relaxation. For the ideal structure, the 
CO/00 state is a very narrow gap insulator. In 
agreement with previous studies, 0, 0| coupling 
to the lattice is however found to be important in 
some respects. Most notably, atomic relaxation at 
the IF shifts the Ti'^^ lower Hubbard band upward 
just enough to lead to conducting behavior, which 
also implies a longer healing length towards bulk be- 
havior, consistent with the experimental indications. 
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